A LEAST SQUARE-TYPE PROCEDURE FOR PARAMETER 
ESTIMATION IN STOCHASTIC DIFFERENTIAL EQUATIONS 
WITH ADDITIVE FRACTIONAL NOISE 



ANDREAS NEUENKIRCH AND SAMY TINDEL 



Abstract. We study a least square- type estimator for an unknown parameter in the 
drift coefficient of a stochastic differential equation with additive fractional noise of Hurst 
parameter H > 1/2. The estimator is based on discrete time observations of the sto- 
chastic differential equation, and using tools from ergodic theory and stochastic analysis 
we derive its strong consistency. 



1. Introduction 

In this article, we will consider the following M'^-valued stochastic differential equation 
(SDK) 

Yt = yo+ h{Y,- §o) ds + T ^ ^ [0^ T]. (1) 

Here uq G M'^ is a given initial condition, B = {B^, . . . , i?™) is an m-dimensional fractional 
Brownian motion (fBm) with Hurst parameter H G (0, 1), the unknown parameter i^Q lies 
in a certain set O which will be specified later on, {b{-]'d), i} G 0} is a family of drift 
coefficients with : M*^ — t- M, and cti, . . . , 0"^ G M*^ are assumed to be known diffusion 

coefficients. 

Let us recall that 5 is a centred Gaussian process defined on a complete probability 
space J-", P). Its law is thus characterized by its covariance function, which is defined 
by 

E{BlBi) = ^-{t'^ + s'^-\t-sn !(,=,), s,teR. 

The variance of the increments of B is then given by 

E\Bi- B'f = \t- s]"^", s,teR, i = l,...,m, 

and this implies that almost surely the fBm paths are 7-Holder continuous for any 7 < if. 
Furthermore, for H = 1/2, fBm coincides with the usual Brownian motion, converting the 
family {B^ , H G (0, 1)} into the most natural generalization of this classical process. In 
the current paper we assume that the Hurst coefficient satisfies H > 1/2 and we focus on 
the estimation of the parameter 'do E Q. Note that the Hurst parameter and the diffusion 
coefficients can be estimated via the quadratic variation of Y, see e.g. [H HI [T5] . 
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Estimators for the unknown parameter in equation ([T]) based on continuous observation 
of Y have been studied e.g. in [21 [HI [T71 [181 EDI |2T1 [22]. Estimators based on discrete time 
data, which are important for practical apphcations, are then obtained via discretization. 
However, to the best of our knowledge no genuine estimators based on discrete time data 
have been analyzed yet. 

We propose here a least square estimator for based on discrete observations of 
the process Y at times {tk] < k < n}. For simplicity, we shall take equally spaced 
observation times with t^+i — = Kn~°' := a„ with given a G (0, 1), k > 0. We call our 
method least square-type procedure, insofar as we consider a quadratic statistics of the 
form 

^ n— 1 

where 5Yu^u2 '■= ~ < Mi < M2 < T and ||a"|p = J2f=i Wjl"^- 

Let us now describe the assumptions under which we shall work, starting from a stan- 
dard hypothesis on the parameter set O: 

Hypothesis 1.1. The set O is compactly embedded in for a given q > 1. 

In order to describe the assumptions on our coefficients b, we will use the following 
notation for partial derivatives: 

Notation 1.2. Let f : M.^ x Q ^ be a C^^'P'^ function forpi,p2 > 1. Then for any tuple 
{ii,. . .ipj e {1, . . ■,dy\ we set di^,„i^J for g^-^^^— ■ Moreover, we will write d^f resp. 

d-af for the J acobi- matrices {d^J, d^J) and {d^J, ...,d^J) . 

With this notation in mind, our drift coefficients and their derivatives will satisfy a 
polynomial growth condition, plus an inward condition which is traditional for estimation 
procedures in the Brownian diffusion case (see e.g [5| [T6]): 

Hypothesis 1.3. We have b E C^'^{W^ x 0;M'^) and there exist constants Ci,C2 > and 
N eN such that: 

(i) For every x,y eM.'^ and d E Q we have 

{b{x; 13) - b{y; -d), x - y) < -Ci\x - y\^ 

(ii) For every x eW^ and E Q the following growth bounds are satisfied: 

<C2(l + |x|^), \dMx;^)\ < C2 (l + , \d^b{x;^)\ < C2 (l + . 

As a consequence of the above assumptions on the drift coefficient and the initial 
condition, for given i^Q E Q the solution of equation ([T]) converges for t — j- 00 to a 
stationary and ergodic stochastic process {Yt,t > 0), see the next section. 

Finally, we also assume that our-drift coefficient is of gradient-type, i.e.: 
Hypothesis 1.4. There exists a function U E C'^'^{M.^ x 0; M) such that 

d^U{x;^) = b{x;^), x G ^9 G 9. 
With those assumptions in mind, we obtain the following convergence result: 
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Theorem 1.5. Assume that the Hypotheses \l.l\ [7751 and are satisfied for equation (J\) 
and that we moreover have H > 1/2. Let Qn{'&) be defined by Then we have 

sup|g„(i9)- (E|6(Fo;^o)l'-E|6(Fo;^)P)| ^0 (3) 

in the P-almost sure sense. 

This convergence is in contrast to the case H = 1/2, i.e. to the case of SDEs with 
additive Brownian noise. There it holds 

sup |g„(^) - (E |6(Fo; ^o) - b{Yo; | ^ (4) 

in the P-almost sure sense, see Remark 13. 6[ and usually the consistent least squares 
estimator 

n-l 

argmin^ge ^ \5Yt^t,+, - 

k=0 

is considered, see e.g. EE]. The difference in the limits and (jl]) is due to the higher 
smoothness and long-range dependence of fractional Brownian motion for H > 1/2. Our 
estimator can thus be seen as a "zero squares" estimator instead of a classical least square 
estimator. In order to show its convergence, we shall work under the following natural 
identifiability assumption: 

Hypothesis 1.6. For any ^ '"^c have 

E|6(Fo;^9o)|' = E|6(Fo;^9)|' ^ff ^ = ^0- 

With this additional Hypothesis, the main result of the current article is the consistency 
of the zero squares estimator based on the statistics 

Theorem 1.7. Assume that the Hypotheses \l.l\ fTT^j \1.4\ and \1.6\ are satisfied for equation 
(Qy and let H > 1/2. Let Qn{'&) be defined by and let dn = argmin^^Q \Qn{'&)\. Then 
for any ^ "we have lim„_!.oo = ^^o in the P-almost sure sense. 

Note that minimizing |(5n(''^)| is of course equivalent to finding the zero of Qn{'&)- Let 
us shortly compare Theorem 11.71 with the existing literature on estimation procedures for 
fBm driven equations: 

(i) Most of the previous results, see e.g. [21 [HI [T71 |2T], deal with the one-dimensional 
fractional Ornstein-Uhlenbeck process in a continuous observation setting. In particular, 
for this process simple continuous time least-square estimators are obtained in |2l[l3], for 
which also covergence rates and asymptotic error distributions are derived. Compared to 
these results our estimation procedure covers a broad class of ergodic multi-dimensional 
equations and relies on discrete data only. 

(ii) A general estimation procedure based on moment matching is established in |2U] . 
However, the main assumption in [2D] is that many independent observations of sample 
paths over a short time interval are available, which is not the case in many practical 
situations where rather one sample path is discretely observed for a long time period. 
Let us also mention the article [3], in which a general discrete data maximum likelihood 
type procedure has been designed for parameter estimation in both the drift and diffusion 
coefficients, however without proof of consistency. 
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(iii) Our current work probably compares best with the maximum hkehhood estimator 
analyzed in [22]. The latter pioneering reference focused on one- dimensional SDEs of the 
form 



with /i : M — 7- R satisfying suitable regularity assumptions. Strong consistency is obtained 
for the continuous time estimator and also for a discretized version of the estimator. 
However, the discretized estimator involves rather complicated operators related to the 
kernel functions arising in the Wiener-integral representation of fBm, which are avoided 
in our approach. Moreover, in contrast to j22] the consistency proof for our estimator 
does not rely on Malliavin calculus methods. 

So, in view of the existing results in the literature. Theorem 11.71 can be seen a step 
towards simple and implementable parameter estimation procedures for SDEs driven by 



Finally, let us comment on the assumptions we have imposed on the drift coefficient 
and on the Hurst parameter: 

(a) The hypothesis of Theorem 11.51 are standard for the case H = 1/2, except Hypothe- 
sis 11.41 which restricts us to gradient-type drift coefficients. We require this condition to 
show an ergodic-type result for weighted sums of the increments of fBm, see Lemma 13.21 
However, this Hypothesis 11.41 is also implicitly present in the additional condition of |16[ 
Theorem 1]. 

(b) It can easily be shown that whenever is a one-dimensional coefficient (namely for q = 
1), Hypothesis 11.61 is satisfied if the drift coefficient is of the form b{x] '&) = §h{x) for some 
/i : M*^ — )■ M'^ and if the stationary solution is non-degenerate, i.e. we have ElFoP 7^ 0. The 
latter conditions hold in particular in the case of the ergodic fractional Ornstein-Uhlenbeck 
process. It would be nice to obtain criteria for richer classes of examples, but this would 
rely on differentiability and non-degeneracy properties of the map d H- 'Ei\h{Y d)^^ (see 
|11| in the Markovian case). We wish to investigate this question in future works. 

(c) Even if the noise enters additively in our equation, we still need the assumption 
H > 1/2 in order to prove Theorem 11.71 Indeed, this hypothesis ensures the conver- 
gence of some deterministic and stochastic Riemann sums in the computations below (see 
Remark 13.31 for further details) . Whether an adaptation of the proposed zero squares 
estimator is also convergent in the case H <l/2 remains an open problem. 

Let us finish this introduction with the simplest example of an equation which satis- 
fies the above assumptions: namely the one-dimensional fractional Ornstein-Uhlenbeck 
process given by 



For t — )■ oo this process converges to the stationary fractional Ornstein-Uhlenbeck process 



dYt = ^QhiYt) dt + dBt 



fBm. 



dYt = ^QYtdt + dBt, 
with -i^o < 0. The solution of this SDE reads as 





t > 0, 
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see e.g. [7|. Here straightforward computations yield the exphcit estimator 



U r. 



E 



n—l -tr2 
k=0 ^tk'^ri 



E 



n-l 
k=0 



\ \ Efc=o 



Otr. 



Notice that even in the case of the fractional Ornstein-Uhlenbeck process, we could neither 
proof the consistency nor show the inconsistency of our estimator for H < 1/2. 

The remainder of this paper is structured as follows: In Section [2] we give some auxiliary 
results on stochastic calculus for fractional Brownian motion. Section |3] is then devoted 
to the proof of our main theorems. 



2. Auxiliary Results 

2.1. Ergodic Properties of the SDE. To deduce the ergodic properties of SDK ([1]) 
we will work without loss of generality on the canonical probability space (^2,^-", P), 
i.e. Q = Co(M, M'") equipped with the compact open topology, J-" is the corresponding 
Borel-(T-algebra and P is the distribution of the fractional Brownian motion B, which is 
consequently given here by the canonical process Bt{u) = uj{t), t G M. Together with the 
shift operators Of : Q ^ Q defined by 

etuj{-) = u{- + t)-u{t), teM, uen, 

the canonical probability space defines an ergodic metric dynamical system, see e.g. [9]. 
In particular, the measure P is invariant to the shift operators 6t, i.e. the shifted process 
{Bs{6f))sm is still an m-dimensional fractional Brownian motion and for any integrable 
random variable F : f2 — )■ M we have 

lim 1 r F{9t{u))dt = EF 

for P-almost all w G f2. Owing to the results in Section 4 of [7] we have the following: 

Theorem 2.1. Let Hypothesis hold. Then for any i^Q E Q the following holds: 

(i) Equation ([1]) admits a unique solution Y in C^(Rj^] Mf^) for all X < H. 

(ii) There exists a random variable Y : Q ^ 'R'^ such that 

lim \Yt{uj)-Y{etUj)\ = 

for P-almost all u E Q. Moreover, we have E\Y\p < oo for all p > 1. 

Note that the law of Y must coincide with the attracting invariant measure for ([T]) 
given in [10], see also [T2| I13|. Moreover, proceeding as in we have: 

Proposition 2.2. Assume Hypothesis \1.3\ holds true. Then for any -i^o G 6 and p > 1 
there exist constants Cp, kp > such that 

E\Yt\P<Cp, E\Yt-Ys\P <kp\t-s\P", for all s,t>0. 

The integrability of Y now implies the ergodicity of equation ([1]): 
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Proposition 2.3. Assume Hypothesis holds true. Then for any i^Q & Q and any 

f e C^(M.'^;R) such that 

+ <c (1 + 1x1^), xeR', 

for some c > 0, N & N, we have 

lim i r f{Y,)dt = Ef{Y) P-a.s. 

Proof. Since the shift operator is ergodic and / has polynomial growth, we have 

lim i r f{Y{9,)) dt = E/(F) P-a.s. 

T^oo T Jq 

Moreover, since 

lim \Yt{uj)-Y{etUj)\ = 

t—^oo 

by Theorem 12.11 and / is polynomially Lipschitz, the assertion easily follows. 

□ 

2.2. Generalized Riemann-Stieltjes Integrals. We set 

||/|U;M= sup |/|a;M= sup 

te[a,b] s,t£[a,b] r ~ ^1 

where / : M ^ M" and A G (0, 1). 

Now, let / e C^{[a, b]; R) and g G C''([a, b]; R) with A + /i > 1. Then it is well known 
that the Riemann-Stieltjes integral j^f{x)dg{x) exists, see e.g. [21]. Also, the classical 
chain rule for the change of variables remains valid, see e.g. [25]: Let / G C^{[a,b]]R) 
with A > 1/2 and F G C^{R] R). Then we have 

F(/(y))-F(/(a))= [' F'{f{x)) df{x), y E [a,b]. (5) 

J a 

Moreover, one has a density type formula: let f,h E C^{[a,b]]R) and g G C'^{[a,b]]R) 
with A + /i > 1. Then for 

ry 

ip:[a,b]^R, ip{y) = j f{x)dg{x), ye[a,b], 

J a 

we have 

h{x) d(p{x) = / h{x)f{x)dg{x). (6) 

J a 

For later use, we also note the following estimate, which can be found e.g. in |24| . 

Proposition 2.4. Let f,g as above. There exists a constant c\^^ (independent of a,b) 
such that 

ph 

< Cx,^,\f\x■[a,b]\9\|^■,[a,b]\b " ^1^"^^ 



{f{s)-f{a))dg{s] 
holds for all a,b E [0, oo). 
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2.3. The Garcia-Rademich-Rumsey Lemma. We will use the following variant of 
the Garcia-Rademich-Rumsey Lemma [8]: 

Lemma 2.5. Let q > 1, a G (1/g, 1) and f : [0, oo) W be a continuous function. Then 
there exists a constant Ca,q > 0, depending only on a, q, such that 



2.4. A Lemma on Pathwise Convergence Rates. The following Lemma (see e.g. 
|19|). which is a direct consequence of the Borel-Cantelli Lemma, allows us to turn con- 
vergence rates in the p-th mean into pathwise convergence rates. 

Lemma 2.6. Let a > 0, po & N and Cp G [0, oo) for p > pq. In addition, let Z^, n G N, 
he a sequence of random variables such that 

(E|Z„ni/^'<Cp-n-" 

for all p > po and all n G N. Then for all e > there exists a random variable rjf, such 
that 

\Zn\ < Ve ■ n-''^' a.s. 
for all n G N. Moreover, 'E\r]^\^ < oo for all p > 1. 



2.5. Quadratic Variations of Fractional Brownian Motion. The following result for 
the behavior of the quadratic variations of a one-dimensional fractional Brownian motion 
P with Hurst parameter H is well known, see e.g [23]. Indeed, for if < 3/4 we have 



^ n— 1 

E|-E0Wi/3r-i] 



k=0 



2 1 

<ch--, 
n 



(7) 



while for if = |, n > 1, it holds 



^ n— 1 

e|-E0W/3P-i] 



k=0 



2 log(ra) 
< C3/4 

n 



Finally if ii G (|, 1) then 



^ n— 1 

E|-E0Wi/3r-i] 



fc=0 



2 1 

< Ch 



n 



(9) 



Here, Ch > ^ denotes a constant depending only on H . 
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3. Proof of Theorems [L5J and [L/7| 

We will denote constants, whose particular value is not important (and which do not 
depend on -d or n) by c, regardless of their value. 

Recall that 

n-l 



For t > 0, setting 

m 

Ft = Y,(^M'\ and n 



|2 11 ||2^2i? 



and moreover using the notation 

S^o-^fK^) = K^'^ ^) - K^'^ ^o), and 5Ft^tk+i = Ft^+i - Ft^, 
it is readily checked that 



fc=0 



fc=0 



^ n— 1 ^ n— 1 



2 



(10) 



n-l 



fc=0 



n.-l 



fc=0 " fc=0 

Note that our assumptions on the drift coefficient imply that 



sup\b{x;'d) -b{y;^)\ < c(l + \x\^ + \y\^) ■ \x - y\ 



for all X, y G M'^ and 



\b{x; A) - b{x; i^a)! < c(l + ■ \A - ^ 



for all X G M and 'di,d2 G O. So, straightforward estimations using Proposition 12.21 give 



Hence for all p > 1 it holds 



and Lemma 12.61 implies 



E 



E 

fc=0 



rk\ 



^ ?^— 1 

lim — 
n^oo na ^ — ^ 



Irfcl = P-a.s. 



k=0 



Using Proposition 12.21 and Lemma 12.61 again, it follows similarly 

n-l 



lim sup — - 



,rk)ar. 



k=0 



P-a.s. 



(12) 
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and, since H > 1/2, we also have 



lim — - 

n^oo naz 



k=0 



P-a.s. 



Plugging relations ( lTT]) -f lT3|) into ( !T0|) . we have obtained that 
where lim„_^oo sup^g0 \Rn{'d)\ = in the P-almost sure sense and 

^ n— 1 _ n— 1 

Til T}/Cy,fi 



(13) 



(14) 



A:=0 



fc=0 



and 



k=0 

The treatment of the terms Qn\'&), Qn\'&) and Qn'' will be carried out in the following 
series of Lemmata. We first show a discrete version of Proposition 12.31 

Lemma 3.1. Let f G C^''^(R.'^ x 6;^^^) be a function such that 

\f{x- ^9)1 < c (1 + , \dj{x- ^)\<c{l + \xf) , \d^f{x- ^9)1 < c (1 + 

for some c> 0, N ^ N, independent o/^? G 0. Then we have 

^ n— 1 

-J2\f{Yt,;n'-WiY;n' 



sup 

i?6e 



fc=0 



P-a.s. 



In particular, we have 



sup |Qi'H^9) - E|5^XF)n ^ P-a.s. 



Proof. Let T„ = na„ and set 



T 



K(^9) = ;^ / \f{Y,;^)\'ds. 



n JO 



The ergodicity of Y yields that there exists a set Ai G J-" with full measure such that 

lim K(^9)(c^)=E|/(F;?9)|2 



for all G 6 n Q'' and all u ^ Ai. The assumptions on / give 



\VnW-VM\<C - [l 

SO Vn is Lipschitz continuous in d and thus 
sup|V;(^?)-E|/(F;^)|2 



T 



n JQ 



(15) 



sup |K(^)-E|/(F; 
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However, from f[T^ and the ergodicity of Y, it also follows that there exists a set A2 E 
with P{A2) = 1 in which the family of random functions : G — > M, ?i G N, is equicon- 
tinuous, and hence the Arzela-Ascoli Theorem yields the desired uniform convergence, 
i.e. 

|2| 



lim sup V;W-E|/(F- 



P-a.s. 



(16) 



Setting 

= \f{Y,-^)\' - \f{Yt,;^)\', t e [h,h+i), k = 0,l, 
it remains to show that 



snp\Gn{t;-t})\dt ^ P-a.s. 

Tn Jo )?Ge 



To this aim, the assumptions on / imply that 



sup \Gn{t; ^)| < c ■ (1 + Yr + • l^t - 



Using Proposition 12.21 and Holder's inequality we obtain 

sup Esup|G„(t;i9)|P < c-<^ 

for all p > 1. Now, Jensen's inequality gives 



(17) 



E 



so ( ]T7j) yields 



^ [ sup \Gn{t; 13)1 dt 

J^n Jo iSGG 



E 



1 Z"^" 
< — / Esup|G„,(t; 
J-n Jo i?ee 



'dt. 



-n Jo i?G© 

for all p > 1 . Lemma 12.61 implies 

Tn 



— I sup \Gn{t]d)\dt 



< c - a. 



pH 



— sup \G nit; ^)\ dt ^ P-a.s. 

Tn Jo i?Ge 



□ 



for n —7- 00. 

We have a similar ergodic result for weighted sums of the increments of the process F. 

Lemma 3.2. Let f G C^'^iW^ x 6; W^) be a function such that 

|/(x; ^)| < c (1 + , \dj{x; ^)\<c{l + \xf) , W{x; ^)\<c{l + \xf) 

for some c > 0, N ^ N, independent of d G O. Assume moreover that there exists a 
function U G C'^^^{W^ x e;M) such that 

d^U{x; ^9) = f{x; ^9), x G ^9 G 6, 

i.e. f is of gradient type. Then, for H > 1/2, we have 

— 5^(/(F*,; ^), 6Ft,t,J + E(6(F; ^0), f{Y; ^)) 



sup 



nar 



k=0 



P-a.s. 



In particular. 



sup\Q^^\iJ) + E{b{Y;^o),5^oAy))\^^ P-a-S. 



1?G0 
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Proof. Let T„ = nan- First note that the chain of variable and density formula for 
Riemann-Stieltjes integrals, see ([5]) and([6]) in Subsection \2.2\ gives that 



1 



T 



Tn 



n Jo 



T 



Tn 



n Jo 



Now the properties of /, Proposition 12.21 and Lemma [2.61 iniplv that 



snp — \U{YT^;d)-Uiyo; 
•See J-n 



P-a.s. 



Moreover, we have 



sup 



T 



n Jo 



P-a.s., 



which can be derived completely analogously to f ll6p . It follows 



sup 



1 



(/(y,;?9),dF„)+E(/(y;^9),6(y;?9o)) 



n JO 

So, it remains to show that 

1 

sup — 



P-a.s. 



^ P-a.s. 



where 



Gr.it- {}) = f{Yt- ^) - f{Yt^- te [4, k = 0,l,.... 



Applying Proposition 12.41 and using the polynomial Lipschitz continuity of / yields, for 
all X<H, 



{Gn{t;^),dFt) 



m n— 1 



A;[tfe;i 



fc+ij ■ 



j=l fc=0 *^[*fc'*fc+il 

From the Garcia- Rademich-Rumsey inequality, see Lemma 12. 5[ and Proposition 12.21 we 
have that 

and also 

Since moreover 



< c - a 



H-X 
n 



1/p 



< c - a 



H-X 



sup (i + iF,r)<c-(F,f 



t£[tfc:tfc + l] 

and supj>QE|yj|^ < oo for all p > 1, it follows that 



( Esup 




f/ 

-'n Jo 


V fee 



{Gn{t;^),dFt) 



< c - a. 



fe+ii 



2H-1 

n 



Now Lemma [2.61 implies f llSp . since H > 1/2. 



(19) 



□ 



Remark 3.3. As mentioned in the introduction, the condition if > 1/2 is used in our 
proofs. Specifically, it is invoked in the convergence of the weighted stochastic integral ([T9|) 
and also in order to derive f[T3l). 
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The following Lemma deals with the remaining term, i.e. the quadratic variations of 
the process F: 



Lemma 3.4. We have 

hm Qf = hm ^ ^ {m,t,J^ - ||af ) = V-a.s. 



fc=0 



with ||cr||2 = Xl^i IcTjP- 
Proof. We have 



tktk+l 



|2 11^112^2// 
I - ll^ll "n 



(j) |2 _ 2/f 



r(2) 



Owing to the scaling property of fBm it follows that 



E 



n-l 
fc=0 



(1) 
k 



a. 



n—l m 



fc=0 j=l 

Since all moments of random variables in a finite Gaussian chaos are equivalent, it follows 
from (ITD-© that 

p 



E 



n-l 



fc=0 



< c - |logrn)P-n^'(2^-^) 



and consequently 



E 



n-l 

E^ 

fc=0 



(1) 



Since a„ = k ■ n with a G (0, 1) Lemma [2.61 now implies that 

n-l 



lim 



n-i.oo nat 



El 

k=0 



(1) 
k 



P-a.s. 



(2) 

So it remains to consider the off-diagonal terms, i.e. . Here we can exploit the following 
trick: Let /3 and /3 be two independent fractional Brownian motions with the same Hurst 
index. From dZ])-® we clearly have that 

n-l 



A;=0 



satisfies 



EWJP <c- al^P ■ I logWP ■ nP^^^^^l 



However, setting B^'^ = (/? + l3)/V2 and B^^) = (/? - l3)/V2, then B^"^ and B^^^ are two 
independent fractional Brownian motions and 

n-l 



k=0 
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Now we can easily conclude that 



lim — - 

n-j-oo nat 



fc=0 



(2) 
k 



P-a.s. 



□ 



Proof of Theorem \1.5[ Let us go back to expression f ITU]) and f lT^ . Applying Lemma [3. II 
and 13.21 we obtain that 



lim sup 

n-)-oo 



(g(i)(^) -2gf (^)) - (E|6(r;^)p-E|6(F;^o)r) 







almost surely. Furthermore, recall that Lemma 13.41 asserts that lim„_!.oo Qrf ^ = almost 
surely. The proof is now finished. 

□ 



Proof of Theorem \1.7[ Let us first recall the following result (pi 116)): 

Proposition 3.5. Assume that the family of random variables Ln{'&), n G N, -i? G O, 
satisfies: 

(1) With probability one, Ln{'d) — )■ Lid) uniformly in d & Q as n oo. 

(2) The limit L is non-random and Li'do) < L{d) for all "i? G 6. 

(3) It holds Li'd) = Li'do) if and only ifd = do- 
Then, we have 

■dn "i^o P-a.s. 

for n ^ oo, where 

Lni'dn) = minL„(?9). 

i?e0 



The strong consistency of the zero squares estimator follows now from Theorem 11.51 and 

□ 



an application of Proposition 13.51 to \Qr 



Remark 3.6. In the case H = 1/2 we have under similar assumptions that 
sup|(g„(i9)-g„(^o)) - (E|6(Fo;i9o)-K>^o;^)|')| ^0 

in the P-almost sure sense, see e.g. Since 

^ n—l ^ n—1 2 n—1 



nat 



k=0 



k=0 



nat 



k=0 



an application of Lemma [3.11 and 13^ (which are also valid for H = 1/2) yield that 



71-1 



lim QniA) = lim — r y^(^^tfctfe+i, '"fc) P-a.s. 
7,— ^oo 71— >oo na •^"^ 



k=0 



However, using the Ito-isometry and the Burkholder-Davis-Gundy inequality we have 



E 



71—1 



k=0 



< C- 



nPa^ 



■Ti+^'/V^/2<c-n-P/2+i-a 
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and Lemma thus gives 

lim Qni^o) = P-a.s. 

Hence we end up with 

sup |g„(^) - (E \b{Yo; ^o) - b{Yo; ^) ^ | -> 0, 

so the hmit of the statistics Qn is different for if = 1/2, where one obtains the standard 
least square estimator. 

References 

A. Begyn (2005): Quadratic variations along irregular subdivisions for Gaussian processes. Electronic 
J. Probab. 10, 691-717. 

R. Belfadli, K. Es-Sebaiy, Y. Ouknine (2011): Parameter Estimation for Fractional Ornstein- 
Uhlenbeck Processes: Non-ergodic Case. Arxiv Preprint. 

A. Chronopoulou, S. Tindel (2011): On inference for fractional differential equations. Arxiv Preprint. 
J.F. CoeurjoUy (2001): Estimating the Parameters of a Fractional Brownian Motion by Discrete 
Variations of its Sample Paths. Stat. Infer. Stock. Process. 4, no. 2, 199-227. 

D. Florens-Zmirou (1989): Approximate discrete-time schemes for statistics of diffusion processes. 
Statistics 20, no. 4, 547-557. 

R. Frydman (1980): A proof of the consistency of maximum likelihood estimators of non-linear 
regression models with autocorellated errors. Econometrica 48, 853-860. 

M. Garrido-Atienza, P. Kloeden, A. Neuenkirch (2009): Discretization of stationary solutions of 
stochastic systems driven by fractional Brownian motion. Appl. Math. Optim. 60, no. 2, 151-172. 
Garcia, A. M., Rodemich, E. and Rumsey Jr., H. (1978). A real variable lemma and the continuity 
of paths of some Gaussian processes. Indiana Math. J. 20, 565-578. 

M. Garrido-Atienza, B. Schmalfuss (2011): Ergodicity of the infinite-dimensional fractional Brownian 
motion. J. Dynam. Differential Equations. 23, no. 3, 671-681. 

M. Hairer (2005): Ergodicity of stochastic differential equations driven by fractional Brownian mo- 
tion. Ann. Probab. 33, no. 2, 703-758. 

M. Hairer, A. Majda (2010): A simple framework to justify linear response theory. N onlinearity 23, 
no. 4, 909-922. 

M. Hairer, A. Ohashi (2007): Ergodicity theory of SDEs with extrinsic memory. Ann. Probab. 35, 
no. 5, 1950-1977. 

M. Hairer, S. Pillai (2011): Ergodicity of hypoelliptic SDEs driven by fractional Brownian motion. 
Ann. Inst. Henri Poincare, Probab. Stat. 47, no. 2, 601-628. 

Y. Hu, D. Nualart (2010): Parameter estimation for fractional Ornstein-Uhlenbeck processes. Stat. 
Prob. Lett. 80, 1030-1038. 

J. Istas, G. Lang (1994): Quadratic variations and estimation of the local Holder index of a Gaussian 
process. Ann. Inst. Poincare 33, 407-436. 

R.A. Kasonga (1988): The consistency of a non-linear least squares estimator from diffusion pro- 
cesses. Stoch. Proc. Appl. 30, 263-275. 

M. Kleptsyna, A. Le Breton (2002): Statistical analysis of the fractional Ornstein-Uhlenbeck type 
process. Stat. Inference Stoch. Process. 5, no. 3, 229-248. 

A. Le Breton (1998): Filtering and parameter estimation in a simple linear system driven by a 
fractional Brownian motion. Stat. Probab. Lett. 38, no. 3, 263-274. 

P. Kloeden, A. Neuenkirch (2007): The pathwise convergence of approximation schemes for stochastic 
differential equations, LMS J. Comp. Math. 10 , 235-253. 

A. Papavasiliou, C. Ladroue (2011): Parameter estimation for rough differential equations. Ann. 
Statist.^ to appear. 

Prakasa Rao, B. L. S. (2010): Statistical inference for fractional diffusion processes. Wiley Series in 
Probability and Statistics, Chichester, John Wiley & Sons. 

C. Tudor, F. Viens (2007): Statistical aspects of the fractional stochastic calculus. Ann. Statist. 35, 
no. 3, 1183-1212. 



LSE FOR SDES WITH ADDITIVE FRACTIONAL NOISE 



15 



[23] C. Tudor, F. Viens (2009): Variations and estimators for self-similarity parameters via Malliavin 

calculus. Ann. Prohah. 37, no. 6, 2093-2134. 
[24] Young, L.C. (1936). An inequality of Holder type connected with Stieltjes integration. Acta Math. 

67 251-282. 

[25[ Zahle, M. (2005). Stochastic differential equations with fractal noise. Math. Nachr. 278, no. 9, 1097- 
1106. 

Andreas Neuenkirch, Fachbereich Mathematik, TU Kaiserslautern, Postfach 3049, D- 
67663 Kaiserslautern, GERMANYneuenkirch@mathematik.uni-kl.de 

Samy Tindel, Institut Elie Cartan Nancy, Universite de Nancy 1, B.P. 239, 54506 
Vandceuvre-les-Nancy Cedex, France tindelSiecn.u-nancy.fr 



